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Abstract. Here we show preliminary calculations of the cooling and contraction of a 2 Mj planet. These calculations, 
which are being extended to 1-10 Afj, differ from other published "cooling tracks" in that they include a core accretion-gas 
capture formation scenario, the leading theory for the formation of gas giant planets. We find that the initial post-accretionary 
intrinsic luminosity of the planet is ~3 times less than previously published models which use arbitrary initial conditions. 
These differences last a few tens of millions of years. Young giant planets are intrinsically fainter than has been previously 
appreciated. We also discuss how uncertainties in atmospheric chemistry and the duration of the formation time of giant 
planets lead to challenges in deriving planetary physical properties from comparison with tabulated model values. 
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1. Introduction 

The direct imaging of extrasolar planets around other stars 
has long been one of the "holy grails" of astronomy. Now, 
ten years after the disovery of the first extrasolar giant planet 
(EGP), this era is apparently now at hand. The two majo r 
discoveries to date are those of IChauvin et^ ( |2004 l2005l) . 
who discovered an object of 5 ± 2 Mj in orbit around a 
25 Mj brown dwarf (2M1207), and lNeuhauser et alJ (l2005l) . 
who discovered a ~ 1 - 42 Mj companion to the T Tauri star 
GQ Lup. Both of these systems are quite young, with ages 
estimated to be -10 Myr for 2M1207 and -2 Myr for GQ 
Lup. 

It must be stressed, however, that these mass estimates 
are based on thermal evolution models, which are most sus- 
pect at young ages. (This is readily apparent from the GQ 
Lup b mass error bar.) Indeed, the most common evolution 
models for the s e brown dwarfs and p l anetary-mass objects 
teurrows etal] 1 19971 IChabrier et all l2000t iBaraffe et all 
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12003 ) use an arbitrary initial condition — the planet is as- 
sumed to be fully convective with an initial arbitrarily large 
radius. This is important because evolution models (or "cool- 
ing tracks") are used to assign masses for these dim objects. 
While giant planets cool quickly and forget their initial con- 
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ditions at ages of billions of years, at ages of —50 Myr or 
less, these initial conditions are important. Using evolution- 
ary models that do not include a valid formation mechanism 
to assign a mass to a young object is quite inaccurate, es- 
pecially at the extremely young ages of these new planetary 
candidates. Modelers have stressed this point. 

This very issue has arisen in t he past few ye ars with the 
study of very young brown dwarfs. iBaraffe et alJ ( 12002 ) have 
analyzed various compilations of evolutionary models of low 
mass stars and brown dwarfs. They also tested various ini- 
tial conditions on their own models and find that the cooling 
tracks differ significantly at young ages (less than —5 Myr), 
which could lead to very different mass or age estimates, 
based on the luminosity of a given object, or group of ob- 
jects. Wuchterl & Tscharnuter (2003) have published forma- 
tion models of stars and brown dwarfs down to 50 A/j. These 
models follow the collapse of an initial cloud all the way to 
the main sequence, while making the minimum number of a 
priori assumptions. These studies highlight the importance of 
initial conditions to the formation of stars and brown dwarfs. 

In tandem with this theoretical work, new observa- 
tions have begun to test the predic tions of evo l ution mod- 
els at young ages. These include ICIose et al ] I2OO5 ) and 
iReiners et alJ ( l2005i) who, with dynamical mass determina- 
tions and age estimates, find that commonly used M dwarf 



2 METHODS 



2. 1 Core Accretion-Gas Capture Model 



evo lutionary tracks (e.g. lBurrows et aL '.'l997^.'Cha brieretall 
I2OOO1) overestimate the l uminosity of these objects at ages of 
^5 - 100 Myr ( for the Close et alJ work) and ^ 5 Myr (for 
the lReiners et alJ workV Mohantv et al. ('2004a'b') have used 
high resolution synthetic spectra to determine Tcs and grav- 
ity for a variety of objects down to perhaps 10 Mj in Up- 
p^r Scorpius, and find disagreement with the Chabrie r et alJ 
11200(1) models. Reiners et al. , in this volume, update the Mo- 
hanty et al. work with new model spectra and find that mod- 
els are consistently overestimating the luminosity at young 
ages. Therefore, as of July 2005, there appears to be dis- 
agreement between these observ ations and m odels (although 
iLuhman et alJ2005l question the lClose et alJ age estimate for 
the object AB Dor). 

To date, there has been little attention paid to the depen- 
dence on formation to the evolutionary history of planets. 
The two leading scenarios for the formation of giant planets 
are th e core accretion-gas capture theory (e.g.. Pollack et al., 
ll£96), which is the generally accepted theory of giant planet 
formation, and the direct gravitational collapse theory (e.g., 
IBossL EoOI: Mayer et al., 2004), which has some proponents. 
Here we investigate the evolution of a 2 AIj giant planet, 
under the assumption that the planet forms via the core 
accretion-gas capture method. We compare the evolutionary 
tracks to those that assume what we term an "initially hot 
start," meaning the planets are assumed to be adiabatic at ev- 
ery age. For these types of models it is common practice to 
pick an initially very hot adiabat, from which point the planet 
cools extremely rapidly, and then plot the cooling history for 
reasonable ages, say from 1 Myr to 10 Gyr. Thes e mode ls are 
typified by [B urrows et al. (Il997l) . iBaraffe et al] (|2003'), and 
iFortnev & Hubbard ( 2004.) . Our models that incorporate the 
core accretion-gas capture formation mechanism are consid- 
erably less luminous at young ages than other published mod- 
els, and these differences can last tens of millions of years, 
not just "a few" million years. 



2. Methods 

2.1. Core Accretion-Gas Capture Model 

The core accretion-gas capture model (a.k.a. core nucle- 
ated accretion) of giant planet formation (Mizuno, 1980; 
iBodenheimer & Pollack[fl986HPollack et allll99dh involves 
the formation of a giant planet by capture of nebula gas by 
a so lid core. These stages are as follows CBodenheimer et al. 



1. Dust particles in the solar nebula form planetesimals 
that accrete into a solid core surrounded by a low mass 
gaseous envelope. Initially, solid runaway accretion oc- 
curs, and the gas accretion rate is much lower than that 
of solids. As the solid material in the feeding zone is de- 
pleted, the solid accretion rate is reduced. 

2. The gas accretion rate steadily increases and eventually 
exceeds the solid accretion rate. The protoplanet contin- 
ues to grow as the gas accretes at a relatively constant 
rate. The mass of the solid core also increases but at a 



slower rate, and eventually the core and envelope masses 
become equal. 

3. Runaway gas accretion occurs and the protoplanet grows 
rapidly. The evolution up to this point is referrred to as 
the nebular stage, because the outer boundary of the pro- 
toplanetary envelope is in contact with the solar nebula 
and the density and temperature at this interface are given 
nebular values. (See Figure^) 

4. The gas accretion rate reaches a limiting value defined by 
the rate at which the nebula can transport gas to the vicin- 
ity of the planet. Subsequent to this point, the equilibrium 
region of the protoplanet contracts inside the effective ac- 
cretion radius, and gas accretes hydrodynamically onto 
this equilibrium region. This part of the evolution is con- 
sidered to be the transition stage. 

5. Accretion is stopped by either the opening of a gap in the 
disk as a consequence of accretion and the tidal effect of 
the planet, or by dissipation of the nebula. Once accretion 
stops, the planet enters the isolation stage. 

6. The planet contracts and cools to the present state at con- 
stant mass. (It is the earliest part of stage 6 that is a focus 
of this article.) 



Simulations of this process have been successful in ex- 
plaining m^n^_o£_the observed characteristics of the giant 
planets jPollack et all 1 199^ . although some important unre- 
solved issues remain. Perhaps chief among these is the abun- 
dance of atmospheric dust grains, which play a major role 
in controlling the times cale of accretion. Recent calculations 
("Hubickyj et al.', '2005^ show that a model in which inter- 
stellar medium (ISM)-like grains are found throughout the 
accreting protoplanetary envelope takes 6 Myr before rapid 
gas accretion onto the planetary core begins. A comparable 
model in which the grain opacity is arbitrarily reduced to 
2% of the ISM value, presumably because of sedimentation 
of grains from high in the atmosphere to their condensation 
level, takes only 2.2 Myr to reach this point. The physical 
mechanism involved is that an atmosphere with grains forms 
a more effective thermal blanket so that the planet retains heat 
from its accretion for a longer time. The planet takes longer to 
contract and thus the onset of rapid gas accretion is delayed. 

Recent calculations bv lPodol^ ll2003h indicate that these 
low opacity cases are likely closer to reality, and the ac- 
tual opacit y may well b e less than the "2% of ISM" used 
by Hubick v! et alJ ( l2005l) . This conclusion is also supported 
by observations of cloudy L-dwarfs, where models that as- 
sume the opacity of grains as given by the ISM distribu- 
tion, without sedimen tation, are a very poor match to the data 
llChabrier et all [2000) because they greatly overestimate the 
dust opacity. How ever, models that allow for g rains to grow 
and sediment (see lAckerman & Marlevi I2OOII) match spec- 
troscopic data very w ell ( Marl ev et all l2002':'Burgasse r et all 
l2002t iMarlev et all .2004) . In short, the timescale argu- 
ment against the core accretion-gas capture scenario is much 
weaker than had been previously thought, when opacities 
more appropriate to protoplanetary atmosphere conditions 
are used. 



2.3 Atmosphere Models 



3 PLANETARY EVOLUTION 





60 

40 
20 




1 




2 






3 


j cool 

\ ^ - 

4 contract — 

] MxY - 
























1 


( (10= yr) 


2 










-2 






















-4 

-6 
-8 




































■ cool 


- 
- 




-10 




1 




2 








& 

contract 















< (10' yr) 


2 










-2 
-4 












U„ 








■S 
\ 

■5 


-6 
-8 
-10 




/ 

f" 1 


Ma 

2 


: cool 
~': contract 













1 


'(10' yr) 


2 










4 






















2 




















\ 






















0:; 

























-2 
-4 




1 




2 








: cool 

contract 





t (10" yr) 



Fig. 1. Ttie mass (units of M®), ttie luminosity (units of L0), 
tlie accretion rates (units of M^/year), and ttie radii (units of 
JJj) ar e plotted as a function of time (units of million years) 
for the iHubickvi et alJ ( l2005h baseline formation model of a 
1 A/j planet. The first three phases of formation of the proto- 
planet (see text) are labeled. Dotted line: phase 1, solid line: 
phase 2, second dotted line: phase 3. The dashed line is con- 
traction and cooling (and is merely illustrative here), at which 
time we switch to an evolution code with fully nongray atmo- 
sphere models. 

2.2. Initial Conditions for Cooling 

Uncertainties in the opacity of the post-formation atmo- 
sphere, when the "isolated" planet begins to cool, lead to un- 
certainties in the luminosity, spectra, and colors of a newly 
formed planet. Although the planet formation code could be 
used to follow the subsequent cooUng of the planets, the 
grain-laden atmospheres that are assumed, while perhaps re- 
alistic during the formation phase when many Earth masses 
of solid planetesimals are ablating into the proto-atmosphere, 
are not correct post-formation. After formation, our proce- 
dure is to switch over to our fully non-gray EGP/brown dwarf 
atmosphere code, so that we can compute a realistic atmo- 
sphere grid for the planets' subsequent evolution. This evolu- 
tion is computed with a sep arate planetary evolution code (see 
iFortnev & Hubbardlbooi . We are also able to use our atmo- 
sphere code to compute spectra and infrared colors, which the 
planetary formation code cannot do. 

Specifically, our procedure is to take the radius (at a pres- 
sure of 1 bar) of the fully formed planet, just after accretion 
has terminated, and use this as the starting point for subse- 



quent evolution. The first model planet calculated with the 
structure/evolution code must match this same radius. (Es- 
sentially no change in the internal struct ure accompanies this 
change in models as both codes use the lSaumon et at Jl 19951 
H/He equation of state.) The planet then contracts and cools 
over time, as the adiabatic interior cools and the atmosphere 
radiates this energy away. These models include a 3.2 g cm~'^ 
solid core which grows during formation and reaches 16 
by the time the planet is fully formed. The exact struc- 
ture of the core has little effect on the evolution, although in 
the future it will be worthwhile to inlude the effects of self- 
compression on the core density. 

2.3. Atmosphere Models 

We have computed a fully nongray grid of model atmo- 
spheres, which serve as the upper boundary condition in 
the evolutionary calculation. This grid spans a range in Toff 
from 85 to 1500 K and in gravity from log g = 2 to 4 
(where g is in cm s^^). To compute our grid we employ 
a model atmosphere code that has been used for a variety 
of planetary and substellar objects. The code was first used 
to generate pressure-temperatur e (P-T) profiles and spectra 
for Titan bv'McKav et al.i ri989l) . It was significantly revised 
to model the atmosphe res of brown dwarfs ( Marlev et aQ, 
a 996; Burrows et al l il997l: iMarlev et all l2002l). U ranus 
dMarlev&McK av. 199<^, coo lEGP s dMarley*, 'l998), and 
"Hot Jupiters" ( Fortnev et al.l l2005b . The basic radiative 
transfer solving scheme was developed bv lToon et alJ ( Il98^ . 

We use the elemental abundance data of Lodders ('2OO3J) 
and compute chemical e quili brium compositions follow ing 
iFedev & LoddersI ill 994 and Lodders & Feelevl ( l2002i) . In 
addition, we maintain a large and constantly updated opacity 
database. Each individual radiative-convective P-T profile is 
arrived at iteratively until the net flux in each layer is con- 
served to at least one part in 10^. Here we neglect the opac- 
ity of silicate and iron clouds, as this study is still somewhat 
provis iona l. The models that w e compare to, burrows et al.l 
(1 19971) and lBaraffeetID(l2003l) . also do not include the opac- 
ity of clouds. 

3. Planetary Evolution 

The main conclusion of this work is shown in Fig. 2. Shown 
is the luminosity of a 2 A/j planet from an age of 10^ 
to 10^ years. The mod els shown are: thic k solid curve 
(this work), dotted curve (Burrows et alJl997 ^. dashed curve 
(Baraffe et at. 2003), and thin solid curve (a variant on our 
work, to be discussed in Section 4.2). Post-formation, at an 
age of 2.5 Myr, our models predict a luminosity only ^^1/3 of 
that predicted by the Burrows et at. and Baraffe et at. hot-start 
models. This difference remains large for tens of millions of 
years, and Burrow s/Baraffe models are still ^^50% more lu- 
minous at an age of 20 Myr. Gradually, the initial conditions 
of the various models are forgotten, and cooling tracks run 
together This is expected, since all of these calculations use 
the same H/He equation of state cSaumon et at- .1995.) and 
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Fig. 2. Luminosity vs. time for a 2 Mj planet for three sets of 
models. The thick solid curve shows our new calculation, in- 
cluding the core accretion-gas capture formation mechanism. 
The planet is ful ly formed at 2 . 2 Myr . The dotted curve is 
for the models of'Burrows et al.' ('1997'). The dashed curve is 
for the models of Baraffe et al. (2003). The thin solid line is 
also our calculation, but the first-post accretion model planet 
is placed at t = 0. The full-width at half-maximum of the 
accretion luminosity (not the log of the luminosity) spike is 
^40,000 years. In reality, the spike is likely to be less sharp 
because of gradual accretion across the gap that the proto- 
planet forms. 

fairly similar atmosphere models. This is a representative cal- 
culation, but our preliminary calculations suggest somewhat 
larger differences at larger masses. However, this will need to 
be investigated further. 

The implications of this calculation for the discovery 
of young Jupiter-mass planets is clear Since hot-start mod- 
els systematically overestimate the luminosity of young gi- 
ant planets, deriving planetary masses based on comparisons 
with cooling tracks, using an object's assumed age and ob- 
served luminosity, will lead one to systematically underesti- 
mate the true masses of the planets. 

4. Discussion 

4.1. Evolution Models 

The implication of our results is that giant planets formed 
by the core accretion-gas capture mechanism are less lumi- 
nous post-accretion than had been previously appreciated. A 
signficant amount of energy is radiated away in the forma- 
tion process. The fully formed planet has a smaller radius at 
young ages than hot start models predict, leading to a lower 
post-formation luminosity. Prior evolutionary model tabula- 
tions overestimate the luminosity of giant planets at young 
ages. The initial conditions for evolution models are not "for- 
gotten" for tens of millions of years. Thus, these planets will 
be more difficult to detect at young ages. 

It is illustrative to compare the integrated radiated energy 
(luminosity x time) of the various models as function of age. 
This is shown in Fig. 3. The line styles are the same as those 



Fig. 3. Integrated radiated energy vs. time for a 2 Mj planet. 
As in Fig. 2, the thick solid curve shows our new calculation, 
including the core accretion-gas capture formati on mecha- 
nism. The dotted curve is for the models of Burrows^etaL 
(|l297f). The dashed curve is for the models of B araffe et alJ 
12003). The dash-dot line is for our own "hot-start" models, 
which are very similar to those are Burrows/Baraffe, but here 
we integrate back to 10 times younger ages (0. 1 Myr). 

in Fig. 2. The location of the lu minosity spike at 2.2 Myr is 
obvious for our calculation. The lBaraffe et alJ P^^ -^^ energy 
curve is always higher than the 'Burrows et al.' {1991) curve 
because as seen in Fig. 2, the luminosity of the Baraffe et al] 
model is larger at every age. The integrated luminosity of our 
new model is greater than either of these models. However, if 
we take our own "hot start" model (dash-dot curve), which 
is very similar to Baraffe/Burrows, and start it at an even 
younger age of 0. 1 Myr, the integrated radiated energy at 1 
Gyr approches that of our core accretion-gas capture model. 
If the starting age for our hot start model is pushed back even 
younger, to 0.01 Myr (not shown), the integrated energy at 
1 Gyr exceeds that of the core accretion-gas capture model. 
(It is unlikely that adiabatic hot start models have relevance 
at these extremely young ages. However, this shows the sen- 
sitiviy of the integrated radiated luminosity to the initial age 
for these models.) The integrated radiative energies are all in 
the range expected, as we find that the gravitational binding 
energy of a 2 Mj planet at 1 Gyr is -8.8 x 10"*^ ergs. 

4.2. Other Detectability Issues 

Another factor that may lead to young Jupiter-mass plan- 
ets being dimmer than previously expected is the result of 
non-equilibrium car bon chemistry. It was first suggested by 
iMarlev et aP ( Il996l) . after the discovery of G1229b, that a 
substantial M-band 4-5 fim flux peak should be a uni- 
versal feature of giant planets and brown dwarfs. In addi- 
tion to the intrinsic emergent flux, this spectral range has 
looked promising for planet detection due to the favorable 
planet/star flux ratio. However, it ha s been known since the 
1970s (see IPrinn & Barshavi 1 19771) fliat the 5 ^m flux in 
Jupiter is less than one predicts from a chemical equilib- 
rium calculation. This is due the dredging up of CO from 
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deeper layers in Jupiter's atmosphere. The chemical con- 
version timescale to CH4 is much slower than the mix- 
ing timescale, ensuring a mixing ratio of CO several or- 
ders of magnitude above the equilibrium mixing ratio. This 
same effect has now been observed in brown dwarf M- 
band photometry ( Golimowski et al., 2004), as anticipated by 
iFeelev & Lodderslj 19961) . This excess CO leads to strong ab- 
sorption at 4.5 /im, leading to diminished flux in M-band 
jSaumon et al., 2003). Consequently, when searching for gi- 
ant planets, L-band may be just as favorable as M-band. 

A factor that seemingly has been ignored to date in the 
detection of young objects around other stars is the actual 
formation time of the planetary candidates. Since the forma- 
tion timescale of giant planets via core accretion-gas capture 
is '^1 - 10 Myr, the assumption of co-evality with the par- 
ent star, which is used for assigning an age to the planet, is 
rather dubious. Interestingly, this effect will work in reverse 
to the intrinsic faintness of the young giant planet — these gi- 
ant planets will be younger than one would assume from co- 
evality — meaning a less massive planet could produce this 
same flux at the given stellar age. It is worth remembering 
that the planetary formation timescale is uncertain, and that 
assuming co-evality at ages of only a few Myr could lead 
to significant errors in determining the physical properties of 
young objects. We n ote that the discr epancy b etween o ur new 
models and those of i Burrows et alJ and Bara ffe et alJ would 
be even greater if our first post-formation planet (here age 2.2 
Myr) was assigned an age of Myr. This is shown as the thin 
solid line in Fig. 2. The luminosity vs. time curve is shifted to 
the left considerably. This further highlights the uncertainties 
in applying evolution models at young ages. 

Finally, we note that given the large physical separa- 
tions of the objects dis covered by£hauvin et al. (2004) and 
iNeuhauser et alJ ll2005l) from their parent stars, it is unlikely 
that these particular objects formed via the core accretion-gas 
capture method. This indeed causes some to claim these ob- 
jects are not "planets" at all, which remains a point of active 
discussion. Whatever the formation mechanism for these ob- 
jects, the point remains that evolution models are uncertain 
at young ages, and that new models and further observational 
calibrations are needed. Our work shows that young Jupiters 
are fainter than has been previously appreciated. In tandem 
with exciting observations, our knowledge about low-mass 
objects at young ages continues to grow. 
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